Theory of topological quantum phase transitions in 3D noncentrosymmetric systems 
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We have constructed a general theory for a new class of topological quantum phase transitions in 3D systems 
with broken inversion symmetry. Although the consideration of the system's codimension predicts the appear- 
ance of a stable metallic phase between the normal and topological insulators, when the band touching occurs 
along the direction with high crystalline symmetry, a direct topological phase transition between two insulators 
can occur. At the quantum critical point (QCP), the energy dispersion becomes quadratic along one direction 
while the dispersions along the other two orthogonal directions are linear, which manifests the zero chirality 
of the band touching point. Due to the anisotropic dispersion at QCP, various thermodynamic and transport 
properties show unusual temperature dependence and anisotropic behaviors. 

PACS numbers: 



The 3D topological insulator (TI) is a new state of mat- 
ter in which the nontrivial topology of bulk electronic wave 
functions guarantees the existence of gapless states on the 
sample's boundary, (J, 0] Because of its topological nature, 
the surface gapless states are protected against small pertur- 
bations preserving the time-reversal symmetry as long as the 
bulk band gap remains finite. Therefore to change the bulk 
topological property, the band gap should be closed at some 
points in the Brillouin zone (BZ) via accidental band cross- 
ings. Recently, such a topological phase transition is realized 
in BiT^Si-zSe^ Ol |4[, by modulating the spin-orbit inter- 
action or the crystal lattice. In the systems with lattice inver- 
sion symmetry such as BiTl(Si_ a: Se 2 ;)2, the topological quan- 
tum phase transition can be described by the (3+1) dimen- 
sional massive Dirac Hamiltonian in general. In this sense, 
the quantum phase transition of 3D TIs provides a new venue 
to study intriguing quantum critical behaviors of 3D particles 
with relativistic dispersion. 15H2I] 

On the other hand, for systems lacking inversion symmetry, 
our understanding of the topological quantum phase transi- 
tion and corresponding quantum critical behavior is still in- 
complete. Based on the consideration of the codimension 
for an accidental band crossing, a stable metallic phase was 
predicted to appear between TI and normal insulator in 3D 
noncentrosymmetric systems. [8] The intervening metallic 
phase between two insulators, dubbed a Weyl semi-metal, has 
topological stability in the sense that there are several gap- 
less points (Weyl points) with nonzero chiral charge at the 
Fermi level. |§] Therefore before every Weyl point is annihi- 
lated by colliding with another Weyl point with opposite chi- 
ral charge, the Weyl semi-metal should stably survive across 
the topological phase transition. In this respect, the recent 
discovery of a direct topological phase transition between 
two insulators in noncentrosymmetric compound BiTel is an 
unexpected surprise. 

EMI At the QCP of BiTel, instead 



of a Weyl semi-metal, several isolated band touching points 
with anisotropic dispersion appear mediating two insulating 
phases, which suggests the diversity of the possible phase dia- 



FIG. 1: (Color online) Generic phase diagrams, resulting from acci- 
dental band touching between conduction and valence bands in the 
3D noncentrosymmetric systems. Three different types of phase dia- 
grams are possible. Here, m indicates an external control parameter. 



grams of noncentrosymmetric systems accessible via acciden- 
tal band crossings. 

In this paper, we propose generic phase diagrams for 3D 
noncentrosymmetric systems that can be achieved through ac- 
cidental band crossings, as depicted in Fig.Q] The key ingredi- 
ent to obtain Fig.QJs the fact that the gap-closing point at QCP 
does not carry any topological quantum number, i.e., the chi- 
rality is zero. Therefore it can be either gapped out leading to 
another insulator (Fig.[T](c)) or split into several Weyl points 
resulting in a Weyl semi-metal. In the latter case, depend- 
ing on whether the trajectory, traversed by the gapless point, 
has a finite torsion or not, the Weyl semi-metal phase persists 
all the way (Fig.Q] (a)) or turns into another insulator (Fig.Q] 
(b)). In all three cases, at the QCP between two neighboring 
phases, the energy dispersion near a band touching point is 
highly anisotropic. Namely, the dispersion is linear along two 
orthogonal directions but is quadratic along the other third di- 
rection. This anisotropic band dispersion induces new power 
laws in the temperature dependence of various measurable 
quantities and anisotropic physical responses. 

Phase transition through accidental band touching.— In 
systems with broken inversion symmetry, the accidental band 
touching between a conduction band and a valence band can 
be described by the following 2x2 Hamiltonian, H (k, m) = 
/o(k, m) + Yli=i fife, m ) T i> where /o,i,2,3 are real functions 
and 71,2,3 are Pauli matrices indicating the two bands. Here 
to describes a parameter, such as external pressure, which can 



be tuned across the phase transition. In particular, we con- 
sider the following situation. For m < m c , the system is fully 
gapped. An isolated band touching point occurs at the critical 
point (k,m) = (k c ,m c ) where A, 2 , 3 (k c ,m c ) = 0. Since f 
does not affect the band touching, we neglect /o in the fol- 
lowing discussion. Then the next question is what happens if 
m is increased further, i.e., when m > m c . To examine the 
system's behavior near the phase transition point at (k c , m c ), 
we derive the low energy Hamiltonian through an expansion 
in powers of q = k — k c and Am = m — m c . Up to the linear 
order of q and Am, f = (/i, /2, fz) T (T stands for transpose) 
can be written as 



f(q, Am) = M q + AmN, 



(1) 



where M,-. 



ait 



If the 



'y — a^"lq=Am=o and Ni = g^| q =A»i 
determinant of M, i.e., DetAI is nonzero, the gap-closing 
condition f = leads to q = — M -1 NAm, which means 
that the gapless point moves as Am varies and persists even 
when Am < contradicting the initial assumption. There- 
fore DetM =0 at the phase transition point. In fact, the sign of 
DetA/ = £ijkMuM2jM3k is the chirality (or chiral charge) 
of the gapless point at Am = 0. Since the chirality is a topo- 
logical number, when a band touching point has a nonzero chi- 
rality, it is stable against small external perturbations. How- 
ever, when DetM =0, it is not topologically protected. There- 
fore when an external parameter changes i.e., Am > 0, the 
band touching point can either be gapped out leading to an- 
other insulating phase or be split into several gapless points 
with zero net chirality generating a stable metallic phase. 
When both of these possibilities are allowed, the insulating 
phase should be preferred since the gapped phase has lower 
energy. 

To understand the nature of the ground state for Am > 0, 
it is useful to rotate the momentum coordinate using a ba- 
sis which manifests the zero chirality of the band touching 
point at Am = 0. Since DetAf=0, M has an eigenvector 
ni with zero eigenvalue satisfying Mm — 0. We introduce 
two additional normalized vectors n 2 , 113, which can form an 
orthonormal basis {ni,na, 113}, and construct a matrix W — 
(ni, n 2 , n 3 ). With the rotated coordinate p = W" _1 q, Eq. (Q~|i 
can be written as, f(p, Am) = U2P2 + U3P3 + AmN, where 
U2,3 = A/n 2 3 . It is worth to note that p\ linear term does not 
appear in f due to the fact that Mm = 0. If the band touching 
occurs at an arbitrary point with no symmetry constraint in the 
momentum space, the leading contribution of pi dependent 
term starts from the quadratic order, which leads to the min- 
imal effective Hamiltonian H(p, Am) = J^ =1 /»(P> Am)r^ 
in which 



f(p, Am) = U2P2 + U3P3 + U4P1 + AmN. 



(2) 



However, if the band touching points are under some symme- 
try constraint, the contributions of other quadratic terms piPj 
should be included to construct the correct low energy Hamil- 
tonian. 





FIG. 2: (Color online) (a) Typical shape of a curve with nonzero 
curvature and torsion in 3D space. 2D Projections of the curve onto 
(b) (t,n) plane, (c) (t,b) plane, (d) (n,b) plane. Here two arrows de- 
scribe the flowing directions of two band touching points which are 
generated from the QCP at the origin. 



Conditions to obtain a semi-metal.— Now let us discuss 
about the ground states for Am > 0. At first we consider 
the case of a stable (topological) metallic phase which can 
be obtained by splitting the band touching point at Am = 
into two gapless points with opposite chiralities. Since there 
are 4 parameters (pi j2;3 and Am) although only 3 condi- 
tions of /1 2,3 = are required to be satisfied to achieve 
a gapless phase, there is a line of gapless points in the (p, 
Am) space in general. Regarding t = Am as a parame- 
ter, the trajectories of the two gapless points form a curve 
p*(i) = (x*(t),y*(t),z*(t)) in 3D momentum space. Then 
the ground state phase diagram can be determined by the 
shape of the curve in the 3D space. Namely, if the trajec- 
tory makes an open curve, the semi-metal phase is expected 
to persist for Am > corresponding to the phase diagram in 
Fig.Q](a). On the other hand, if the curve forms a closed loop, 
there is another transition from the semi-metal to an insulator 
leading to the phase diagram in Fig.Q](b). 

A convenient way to describe a curve C(t) = 
(x(t),y(t), z{t)) is to use an orthogonal coordinate system 
(t, n, b) where t, n, b = t x n are the tangential, normal, 
and binormal vectors, respectively. The shape of a curve C 
near the initial point at t — can be determined once the cur- 
vature kc and the torsion tq are specified at t = 0. II 1311 Here 
kc measures the rate at which the tangential vector changes 
and tc measures the rate at which the curve C deviates from 
being a planar curve lying on the (t, n) space. The relation 
between the shape of a curve and its curvature and torsion can 
be understood by considering two dimensional projections of 
the curve on the planes in which two vectors among {t, n, b} 
are adopted as a basis. After choosing a coordinate properly, 
the projections of the curve C are described by y = ^&x 2 , 



and z 2 



— 2 T c „,3 



y , which are lying on (t, n) 



(t,b), and (n,b) planes, respectively. A typical example of 
a curve in 3D and its 2D projections are shown in Fig. [2] 
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Considering the projections of the curve on (t,b) and (n,b) 
planes, shown in Fig. [2] (c) and (d), respectively, we can see 
that the trajectory cannot form a closed curve under smooth 
variations as long as the torsion tq is nonzero. The semi- 
metallic phase should persist in this case. (Fig.Q](a)) An ex- 
ample where gapless points form a curve with nonzero tc is 
shown in the supplementary material. 11411 On the other hand 
when tc = 0, the curve moves only on the 2D space spanned 
by (t, n). The trajectory of the curve on (t, n) plane shown in 
Fig.|2](&) can make a closed loop as long as the curvature kc 
is finite. In this case, it is expected that there can be another 
phase transition to an insulator through pair annihilations of 
band touching points. (Fig. Q] (b)) Therefore the fate of the 
semi-metallic phase is determined by the torsion tc of the 
curve corresponding to the trajectory of band touching points. 

Conditions to obtain an insulator. - Now we turn to the con- 
ditions to have a gapped phase for Am > corresponding to 
Fig. [2(c). Since the system is gapped for Am ^ 0, the con- 
duction (valence) band should have a well-defined dispersion 
minimum (maximum) near p = 0, which can allow a band 
touching locally in momentum space at Am = 0. Consid- 
ering the minimal 2x2 Hamiltonian with /y^^Am) in 
Eq. (fJJ, the energy eigenvalue of the conduction band is given 

by£' c (p, Am) = ^Jj2i=i 2 3 fi(P> Am). The condition to 
have an extremum for small Am 7^ leads to the following 
three equations 9.; = dE c (p, Am) /dpi = (i = 1,2,3). 
After solving the 3 coupled equations, the general solution 
for the location of the dispersion minimum is obtained as 
pinm _ ^ A/± m ^ BAm) where A, B are some constants. It 
is important to notice that the conduction (valence) band min- 
imum (maximum) should move along the line that satisfies 
pi = and pi — jjPz for both Am < and Am > 
across the phase transition, which can be achievable generally 
when the system has high crystalline symmetry along the line. 
Moreover, due to time reversal invariance, the band structure 
at momenta ±k should be equivalent. Therefore the extrema 
in band dispersions of the gapped phases should move along 
the line with high crystalline symmetry which passes through 
a time reversal invariant momentum at the center. 

To understand the shape of the band near the dispersion 
minimum, we compute the Hessian matrix -ff,™ 11 = g p .^. , 

which has a block diagonal form with H™" = ^13" = at 
p = p min . The other nonzero components of H min satisfies 



#57 = c u Am, Det 



H min 
22 

tt mm 

H 32 



Hmin 
23 

^33 



>o, 



where en is a constant. Interestingly, H™ n changes the sign 
across the phase transition because it is linearly proportional 
to Am. For c\\ < (en > 0), the conduction band has 
a dispersion minimum in all three directions for Am < 
(Am > 0) while it has a saddle point with a negative curvature 
along the p\ direction for Am > (Am < 0). Therefore 
when there is a direct phase transition between two insulators, 
one of them has to possess a saddle point at the bottom (top) 



of the conduction (valence) band along the p\ direction where 
the energy dispersion is quadratic at QCP. 

Topological phase transition— A direct phase transition be- 
tween two insulators via accidental band touching can accom- 
pany the change of bulk topological properties. In 3D time- 
reversal invariant systems, band insulators can be classified by 
Z2 topological numbers Vi (i = 0, 1, 2, 3). HI [l5i [l6ll In the 
BZ, there are 3 pairs of parallel planes, in which k • = 
or tt. (i = 1,2,3) Here 81,2,3 are primitive lattice vectors. 
Since each of the planes satisfies the time reversal symmetry 
Tif(k)T _1 = H(-k) with time reversal operator T, a 2D 
Z2 invariant a? (af ) can be assigned to each plane satisfy- 
ing k • = (k • a t = tt). Since ot\ + a\ = a% + = 
a§ + , only four 2D invariants are independent and deter- 
mine the Z2 invariants of the 3D system in the following way, 
(i/q, v\, v%, 1/3) = {a\-\- a\ , a™, a%, ctg). The strong invariant 
vq distinguishes a TI (i/q = 1) and a band insulator (vq = 0). 
Since vo = a® + otj for any i — 1, 2, 3, if one of 2D Z2 in- 
variants changes by 1 through the accidental band touching, a 
topological phase transition occurs. 

In a 2D BZ with time-reversal invariance, the Z2 invariant a 
is given by the Chern number (modulo 2), which is the integral 
of the Berry curvature over the half of the BZ (with additional 
contraction procedures). B15I1 Therefore if the band touching 
between the valence and conduction bands, which changes the 
Chern number of each band by ±1 per single touching ifnll . 
occurs odd number of times in a half of the BZ, a changes 
by 1 leading to the 3D topological phase transition. There- 
fore when the direct phase transition between two insulators 
happens, if the line with high crystalline symmetry embracing 
QCPs is lying on the 2D planes with time-reversal invariance 
and the number of such lines is odd, the topological property 
of the two insulators should be different. 

Application to the topological quantum phase transition of 
pressured BiTel.- Recently, it is theoretically demonstrated 
that when an external pressure P is applied to BiTel, a direct 
phase transition from a normal insulator to a TI occurs medi- 
ated by accidental band touching points at the critical pressure 
P c 1 1 211 - The origin of the direct phase transition in BiTel can 
be understood based on the general considerations developed 
in previous discussions. The key ingredient is the fact that the 
band touching occurs along the line with high symmetry (A- 
H line) in k z = tt plane (BZ of BiTel is shown in Fig. 1 of 
1 1210. Because of the point group symmetry, the conduc- 
tion (valence) band with the Rashba-type spin-orbit coupling 
develops a dispersion minimum (maximum) along this line 
for any pressure across the accidental band crossing, which 
satisfies the necessary condition to achieve the direct transi- 
tion. Moreover, since the high symmetry line is on the k z — tt 
plane satisfying the time-reversal invariance, the band touch- 
ing can induce the topological phase transition. In Fig. [3] we 
plot the evolution of the band dispersion across the topologi- 
cal phase transition near one of band touching points using the 
band structure obtained by first principle calculations. Ill 811 At 
the QCR the band dispersion is quadratic along one direction 
(pi) and linear along the other two orthogonal directions (p2 
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FIG. 3: (Color online) Evolution of the band structure, obtained from 
first-principle calculations, across the topological phase transition in 
BiTel. Energy dispersion of the conduction/valence bands near one 
of the band touching point in (pi,p2) plane, which is normal to the 
high symmetry line embracing QCPs, is shown for (a) P < P c , (b) 
P = P c , (c) P > P c , respectively. Energy dispersions along the pi 
(red), p2 (green), p$ (blue) directions are shown for (d) P < P c , (e) 
P — P c , and (f) P > P c , respectively. 





D(e) 


C V (T) 




K d(T) 


(TDC(T) 


Weyl semi-metal 


e 1 






InT 


T 


At the QCP 


£ 3/2 




rp3/2 


T -l/2 


jiX/2 



TABLE I: Summary of the leading temperature (or energy) depen- 
dence of various physical quantities for a 3D Weyl semi-metal and at 
the QCP. D(e), Cv, K, xd, <jdc are density of states, specific heat, 
compressibility, diamagnetic susceptibility, and the DC conductivity, 
respectively. For both the 3D Weyl semi-metal and at the QCP, we set 
the chemical potential to be at the band touching point. To compute 
odc(T), we have used the T-linear scattering rate due to Coulomb 
interaction, which is suitable for quantum critical systems. 



and 7)3). Moveover, across the QCP, the band dispersion of the 
insulating phase for P > P c possesses a saddle point, which 
proves the direct topological phase transition. 

Thermodynamic properties at QCP.- The most important 
characteristics of the QCP is the anisotropic dispersion, which 
leads to the following minimal Hamiltonian at the QCP, 



-ffocp(p) = Ap\r\ + vp 2 T 2 + vp 3 T 3 . 



(3) 



where v is the velocity and A is the inverse mass along the 
pi direction. The anisotropic dispersion gives rise to a pecu- 
liar power law in the density of states such as D{e) oc e 3 / 2 , 
which is quite distinct from that for a 3D Weyl semi-metal 
(D(e) oc e 2 ) or a 3D normal metal with quadratic dis- 
persion (D(e) oc e 1 / 2 ). This distinct power law of D(e) 
directly leads to new exponents in the temperature depen- 
dence of various thermodynamic quantities such as the spe- 
cific heat (Cv) and compressibility (k) as summarized in Ta- 
ble U The diamagnetic susceptibility \d also shows an un- 
expected singular behavior. We have computed \d using 
the Fukuyama formula for the orbital susceptibility xd = 



^vT.n^\GlaG lb G la G lb ]. H Here 7a = ^gf and 
a, b are two orthogonal directions perpendicular to the applied 
magnetic field. For the Hamiltonian in Eq. (01, \d is given by 
Xd(0) = cos 2 6>xi + sin 2 9x2, in which xi ~ CiT~ 1/2 and 
X2 ~ x^+^T 1 ^ 2 with C1.2 constants. Here 8 is the angle 
between the external magnetic field and p\ direction. There- 
fore xd shows unusual singular temperature dependence in 
low temperature given by xd ~ T -1 / 2 irrespective of mag- 
netic field directions. 

Anisotropic DC conductivity.- The anisotropic dispersion at 
QCP also induces anisotropic temperature dependence of the 
DC conductivities. Assuming momentum independence of 
the scattering rate t(oj), a straightforward calculation of the 
conductivity tensor using Kubo formula gives rise to the fol- 
lowing expression of the DC conductivities, 



*n(T) = 

022,33 = 



2e 2 \/A 

7lT 2 V 2 

9e 2 



(lid 



M 5/2 (- 



20ir 2 VA 



du\u) 



3/2 



doj) 



(4) 



When the Coulomb scattering is the dominant scattering 
mechanism, we can take - = a 2 T with a = — the typ- 
ical scattering rate of quantum critical systems, taking into 
account of the fact that the low temperature transport is domi- 
nated by the linear dispersion. In this case, the temperature de- 
pendence of the DC conductivity is given by ctu(T) oc T 3 / 2 
and (J22.33 (T) oc T 1 / 2 . On the other hand, when the scat- 
tering due to random potentials dominates the transport, us- 
ing Born approximation, the leading contribution to the scat- 
tering rate can be obtained by rj 2tt^oD(w) with 

7o = "' 2 ° ■ 12011 Here Vb is the impurity scattering poten- 
tial, rii is the impurity density. Then using Eq. (O, we obtain 

9e 2 u 2 „ (rp\ _ 2e 2 A, 



CT33(T) 



207T7Q 



a n (T) = ff^(21n2)T, which also shows 



the anisotropic T dependence. 112111 In fact, Eq. (0} implies 
that, as long as the scattering rate does not show any momen- 
tum dependence, 033 (T)/a%x (T) shows unusual temperature 
dependence, independent of the scattering mechanism such as 

£ag = G,£Twhere<7 «1.8. 

Stability of QCP.- Finally, let us discuss about the stabil- 
ity of the anisotropic QCP in Eq. ([3]) under the influence of 
disorder and long-ranged Coulomb interaction. The effective 
action of the QCP including both random disorder potential 
and 1/r Coulomb interaction can be written as 

S= J dtd 3 x{ij f (idt + Ad 2 T 1 + ^ ivd^tp + F^M^] 



J=2,3 
2 



+ / dtd 3 xd 3 x'{^^) x , t 9 _ x „ (VW.i- 



(5) 



where Vi(x) is a random potential coupled to fermion field 
tp(x) via a matrix Mi. g 2 = e 2 /e where e and e are the elec- 
tric charge and dielectric constant, respectively. We take a ran- 
dom disorder potential with Gaussian invariance whose impu- 
rity average satisfies (T^(x)Vj(x')) = Ay<5^(x — x'). The 
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key characteristics of the Gaussian fixed point in Eq. (O is the 
invariance of the Hamiltonian under the anisotropic scaling of 
spatial coordinates, i.e., X\ — Xi/b 1 / 2 , x 2j 3 = ^2,3/^ accom- 
panied by t = t/b where the tilde indicates the new scaled 
coordinates. Under this scale transformation, A,j transforms 
as Aij = 6~ 1//2 Ajj showing that the disorder is an irrele- 
vant perturbation. Similarly, it is straightforward to show that 
g 2 = b~ 1/>2 g 2 . Therefore both the random disorder poten- 
tial and Coulomb interaction are irrelevant perturbations at 
the Gaussian fixed point. The same conclusion can be ob- 
tained even when we consider a generalized scaling transfor- 
mation, i.e., t = t/b z , X\ — xi/b Zl , and x 2 ,3 — x 2,z/b within 
one-loop perturbative renormalization group calculation. lfl4ll 
Therefore the unusual power laws in various thermodynamic 
and transport properties, which are predicted based on the free 
particle Hamiltonian in Eq. (01 should persist even under the 
influence of the disorder and Coulomb interaction. 
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SUPPLEMENTARY MATERIAL 
1. A curve in 3D and its curvature and torsion 

In this section, we describe the relation between the shape 
of a curve in 3D space and its curvature and torsion. Here we 
basically follow the contents in Ref. [TJ. A convenient way 
to describe a curve C = C(t) — (x(t),y(t), z(t)) is to use 
an orthogonal coordinate system (t, n, b) where t, n, b are 
the tangential, normal, and binormal vectors, respectively. To 
define t, n, b we first consider the arclength, which is defined 
as 



s(i) = / dt' 



dC{t') 



df 



(7) 



For the following discussion, we reparametrize the curve C 
using the arclength s, which makes C to be a function of s, 
i.e.,C = C(s). Then t, n, b are given by 



dC 
ds 



d 2 C , 



d 2 C 



ds 2 



b = t x n, 



and the curvature kq and torsion re are defined as 



kc(s) = 



db 

- = -rc(s)n. 



(8) 



(9) 



Therefore nc measures the rate at which the tangential vector 
changes and tc measures the rate at which the curve C devi- 
ates from being a planar curve lying on the (t, n) space. In the 
case of the torsion tq, it has a following equivalent expres- 



sion, 



TC 



J_/dC (fC_ d 3 C \ 
K 2 \ds X ds 2 ' ds 3 /' 



(10) 



where (A, B) indicates the inner product of two vectors A and 
B. 

The relation between the shape of a curve and its curvature 
and torsion can be understood by considering two dimensional 
projections of the curve on the planes in which two vectors 
among {t, n, b} are adopted as a basis. To represent an arbi- 
trary curve C(s) we can choose a coordinate which satisfies 

C(0) 



= 0, m 

as 



s=0 



(1,0,0), and 



Then from Eq. (HOI , we obtain 
where a% t 2 are arbitrary numbers 



ds 3 



s=0 



s=0 



= (0,k o ,0). 
(ai 7 a 2 ,KcTc) 
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Now we consider the Taylor expansion of C(s) 
(x(s),y(s),z{s)) at s = 0, 



C(s) = C(0) + 



(is 



s 2 d 2 C 




s 3 d 3 C 


o + 2 ds 2 


H 




6 g?s 3 



+ 0(s 4 ), 

(11) 



which gives rise to 

x(s) = s + 0(s 3 ), 

,( S ) = ^ S 3 + 0( S 4 ) . 



(12) 



From the leading order terms, the projections of the curve C 

are described by y = ^f-x 2 , z — Kc ^ c a: 3 , and z 2 — f^-y 3 , 
which are lying on (t,n), (t,b), and (n,b) planes, respec- 
tively. A typical example of a curve in 3D and its 2D pro- 
jections are shown in Fig. [2] It is to be noted that the projec- 
tions of the curve on (t, b) and (n, b) planes, shown in Fig. |2] 
(c) and (d), respectively, cannot form a closed curve under 
smooth variations of the curve. Since the shapes of these 2D 
projections remain the same as long as the torsion is nonzero 
at the origin, we obtain the zero torsion condition tq = to 
form a closed curve. On the other band, when tq = 0, the 
curve moves on the 2D space spanned by (t,n). The trajec- 
tory of the curve on (t, n) plane shown in Fig. |2] (£») can make 
a closed loop as long as the curvature kc is finite. 



2. Effective Hamiltonian for topological phase transition in 
BiTel 

Let us first consider a general Hamiltonian defined in mo- 
mentum space H(k) = H(k x , k y , k z ). Then the combined 
symmetry operation TM, which is the combination of the 
time reversal (T) and mirror (M : y — > —y), imposes the 
following constraint to the H(k), 

{TM)H(k x ,k y ,k z )(TM)- 1 
= H {k x , ky, kz) = H( k x , ky, k z ). (13) 

Assuming that the band touching point exists at k c = 
(0, ky^ c , 7r), we derive the low energy Hamiltonian consider- 
ing small momentum deviation from the band touching point. 
Due to the TM symmetry, the effective 2x2 Hamiltonian 
-H2x2(q) with q = k k c , satisfies the following constraint, 

(TM)H 2X 2(qx,q v ,qz)(TMy 1 = H% x 2 (q X) q v , &) 

= ^2x2(^11 ky ky,c, k z — 71") 

= H 2X 2(-k x , k y - k VtC , -k z - 7r) 
= H 2X 2(-k x , k y - ky, c , ~k z + 7r) 

= H 2x2 (-q x ,q y ,-q z ). (14) 

Therefore 



For the generic 2x2 Hamiltonian H 2x2 (q) given by, 

H 2x2 (q) = /i(q)ri + / 2 (q)r 2 + / 3 (q)r 3 , (16) 

the constraint in Eq. (15[ leads to the following constraints to 

/i,2,3(q), 

fi,a(-Qx,Q y ,-qg) = fx, z (qx,q v ,qz) 
fa(-q x , q y , -q z ) = -f y (q x , q v , q z ) (17) 

which means that /i, 3 (f 2 ) are even (odd) under the simulta- 
neous sign change of q x and q z . Now we expand 71,2,3 near 
the band touching point at q = and P = P c in the powers of 
q x ,y,z an d = P — Pc- Due to the symmetry constraint in 
Eq. (fTTI i. up to linear order in q x ,y,z an( l = P ~ Pc, 1 1,2,3 
are given by 



fa =N 1 AP + M 12 q y 
fa = M 21 q x + M 23 q z 
h = N 3 AP + M 32 q y 



(18) 



where jVi,3 and M12, 21,23,32 are constants. The Hamiltonian 
at the critical point AP = can be written as 

#2x2(q) = M 12 q y Ti + (M 21 q x + M 23 q z )r 2 + M 32 q y T 3 , 

(19) 

Since DetA/ = 0, M always has an eigenvector £1 with 

= (-M 23 ,0,M 21 ) 



zero eigenvalue. Explicitly, £\ — , 

V M 21+ M 23 

where the superscript t means the transpose of a vector. 



Let us introduce £ 2 



= (M 21 ,0,M 23 ) and£* 



(0, 1,0). Then £1,2,3 constitute an orthogonal basis. Using 
this, we consider following linear transformation, 




(20) 



Applying the above linear transformation, the Hamiltonian is 
given by H 2x2 (p) = X«=i,2,3 />(P) T » in which 

h = N X AP + M 12 p 3 , 



fa = V M 2 2 i + M 2 3 p 2 , 
fa = N 3 AP + M 32 p 3 . 



(21) 



It is to be noticed that p\ does not appear in the Hamiltonian 
due to the fact the £1 is the eigenvector of M with zero eigen- 
value. 

To fully account for the phase transition, the terms 
quadratic in pi are necessary. In terms of the rotated momen- 
tum p, the symmetry constraint in Eq. (TTTb can be written as 



H 2x2 {qx,q y ,qz) = H 2x2 (~q x ,q y ,-q z ) (15) 



H 2 X2 (P1,P2,P3) = H 2x2 (-px, ~p 2 ,p 3 ). 



(22) 




(a) 



(b) 



Pi 




(c) 



FIG. 4: (Color online) The hexagonal Brillouin zone of BiTel. 



Collecting terms satisfying the constraint above up to the 
quadratic order in pi, f\ 2, 3 can be written as 

fx = TViAP + M12P3 + aip\ + a 2 p\ + a z p\ + a 4 piP2, 



h = y M 21 + M 23P2 + b 5P2P3 + hP3Pl, 

f 3 = N3AP + M32P3 + Cip\ + C2P2 + C3P3 + c 4PlP2, 

(23) 

where a^, bi, a (i — 1, 2, 6) are constants. 



3. Topological phase transition in BiTel 

Recently, it is shown that when external pressure is applied 
to BiTel, a direct phase transition from a normal insulator to a 
TI occurs mediated by accidental band touching points at the 
critical pressure P = P c . The nature of the topological phase 
transition in BiTel can be understood based on the general 
considerations developed in previous discussions. The key in- 
gredient to understand the phase transition in this material is 
the fact that in the BZ shown in Fig. [4] the band touching oc- 
curs along the A — H direction in k z = tt plane, along which 
the Hamiltonian P(k) is invariant under the combination of 
the time reversal (T) and mirror (M) symmetries. Picking 
one of the A — H direction along the k y axis, the combined 
symmetry operation TM imposes the following constraint 
to the Hamiltonian, H*(k x , k y , k z ) = H(—k x , k y , — k z ) be- 
cause the mirror M changes k y to — k y . This symmetry con- 
straint restricts the structure of the low energy Hamiltonian 
near the gap-closing points. Explicitly, the effective Hamilto- 
nian near one of the band touching points can be written as 
#(p, AP) = £i=i fa, AP) n in which 

/i = NxAP + M 12 p3 + aipl + a 2 p\ + a, z p\ + a 4 piP2, 



h = y M 21 + ^23^2 + hp 2 p3 + &6P3P1, 
h = N 3 AP + M 3 2P3 + Cip\ + C2p\ 



C3V\ 



C4P1P2, 

(24) 



where a,i, bi, a (i — 1, 2, 6) are constants and AP = 
P — P c . p is the rotated momentum coordinates adapted to 
manifest the zero chirality of the band touching point. There- 
fore the terms linear in p\ do not appear in /i,2,3- The detailed 



FIG. 5: (Color online) The 2D projections of the curve satisfying the 
gap-closing condition for the trajectory in Eq. ( 125b for a > 0, /3 > 0, 

7 > 0. 



procedures to derive the above Hamiltonian is shown in the 
Supplementary Material. 

Let us first check the conditions to obtain a semi-metallic 
phase by finding the solution of /i 2,3 = 0. When |p| <C 1, 
h = leads to, P2 oc P3Pi- Inserting this result to the con- 
ditions /1.2 = 0, the general solution of the gapless point has 
the following structure, p\ = ±oVAP, p* 2 = ±^(AP) 3/2 , 
P3 = 7AP with a, f3, 7 constants. Here we assume that the 
gapless point exists only for AP > 0. Therefore the trajectory 
of the gapless points is given by 



P 



/ AP,±/3(AP) 3 / 2 , 7 AP) 



whose curvature and torsion at AP = is given by kq 



(25) 



_ 27 



and tc = -§^-, respectively. Since the curve has a finite tor- 
sion, we expect the trajectory of the gapless points would not 
form a closed loop under the smooth variation of the system, 
which is supported by the corresponding 2D projections of the 
curve shown in Fig. [5] Therefore once a semi-metallic phase 
occurs by splitting the band touching point at AP = 0, the 
gapless point should persist for all AP, which is not consis- 
tent with the predictions of the first principle calculation. 

The only way to describe the transition between two insu- 
lators is the occurrence of the direct transition between them 
via band touching at a single critical point. From the condi- 
tion that the conduction band has a minimum, we have three 



equations of 



dE c (p) 



(i = 1, 2, 3). After careful examina- 



tion of the three coupled equations, we have found that there 
is a unique solution given by 



= (0,0,- 



M\ 2 



AP). 



(26) 



Therefore the direct transition between two gapped phases re- 
quires that the extrema of the valence and conduction bands 
move along a particular direction, the P3 direction through- 
out the phase transition. Interestingly, the P3 direction corre- 
sponds to the k y direction of the original coordinate, which is 
nothing but one of the A — H direction in the BZ. Because 
of the C3 V point group symmetry, the conduction (valence) 
band with the Rashba-type spin-orbit coupling develops a dis- 
persion minimum (maximum) along the A — H line for any 
AP ^ 0. Moreover, since the A - H line is on the k z — n 
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plane satisfying the time-reversal invariance, the band touch- 
ing can induce the topological phase transition. Explicitly, 
in the k z = it plane, the effective Hamiltonian near a gap- 
closing point can be written as a two dimensional massive 
Dirac Hamiltonian, in which AP plays the role of the mass 
term. Therefore a band touching point can change the Z2 in- 
variant on the k z = tt plane by 1 through the sign change of 
AP. Since there are 3 pairs of band touching points on the 
same plane, the strong index Uq changes by 1 reflecting the 
emergence of the TI through the band touching. 

has a block diagonal 



The Hessian matrix H 



d 2 E a 
dpidpj 



form with jfjf = H% 3 m = at p 
components of H mm satisfies 



rj nun 

^33 



> 0, Det 



rrmin 

H n 

rrmin 
H 21 



rrmin 
H 12 
rrmin 
H 22 



The other nonzero 



= CAP, 



where C is a constant. Therefore the conduction band always 
has a minimum along the p 3 direction. On the other hand, the 
determinant of the Hessian matrix for the p\, P2 directions, 
normal to the A — H line, shows sign change across the phase 
transition, which predicts the appearance of a saddle point in 
{pi,P2) plane for one of the gapped phase. Since the normal 
insulator for AP < possesses the conventional dispersion 
minimum or maximum, we can set C < 0. Therefore the TI 
should possess saddle points in energy dispersion. The change 
of the band dispersion across the topological phase transition 
is explicitly shown in Fig.[3]from the band structure obtained 
by first principle calculation. 

4. Perturbative renormalization group approach to the 
quantum critical point 

The Hamiltonian for the anisotropic QCP can be written as 

H = —ivd\T\ — iv&2T2 — Ad 3 r 3 , (27) 

Compared to the Hamiltonian in the main text, the coordinates 
x\ and x 3 are interchanged for convenience. In the main text, 
we have considered the following scale transformations, 



t = bt, x\ 



■I'3 



b l ' 2 x 



x 3 , 



(28) 



which lead to the conclusion that both the point disorder and 
1 jr Coulomb potential are irrelevant perturbation in tree level 
analysis. Here we consider a more general scale transforma- 
tion without specifying the scaling dimension of x 3 and t. 
Namely, 



t = b Z t, X\2 = bX\2t X 3 — b Z3 X 3 . 



(29) 



The effective action for free electrons at QCP can be written 

as 

S = J dtd 3 x[^{id t + ivdm + ivd 2 r 2 + Ad$T 3 )i/>]. 



(30) 



The invariance of So under the scale transformation in 
Eq. (|29l leads to the following scaling relations at tree level, 



4>(x,t) = b 1+ ^ijj(x,t), v 



A 



A. 
(31) 



If we assume that the velocity v is a marginal coupling at tree 
level, we obtain z = 1. In this case, from the scale trans- 
formation of A, we can impose z — 2z 3 < 0, which leads 
to Z3 > 1/2. On the other hand, if we assume that A is a 
marginal coupling, we obtain z = 2z%. From the scaling di- 
mension of v, i.e., z — 1 = 2^3 — 1 < 0, we obtain z 3 < 1/2. 
Therefore z and z 3 should satisfy either z — 1 and Z3 > 1/2 
or z = 2z 3 and < z 3 < 1/2. 

Under this scale transformation, the disorder potential for 
the short-ranged random disorder A^ transforms as A^ = 
b 2z ~ 2 ~ Za Aij, in which the scaling dimension satisfies 2z — 
2 — z 3 < —1/2. Namely, the disorder is always irrelevant. 
Therefore in the following we neglect the disorder effect and 
focus on the influence of 1/r Coulomb potential to the Gaus- 
sian fixed point. 

To perform the scaling analysis and one-loop renormaliza- 
tion group calculation, it is convenient to write the effective 
action in the momentum space, which is given by 



S — Sq + Sc, 



(32) 



in which 



V 



(2tt) 3 



dtd 3 k^[id t - u(fc]Ti + k 2 T 2 ) - Aklr^ip, 
1 



„2 y2 

Sc = ^ — ~ — 777 / dtd 3 qd 3 k 1 d 3 k 2 



9x + ll 



2 (2tt) 9 

x V f (ki +q,^(k 1 ,t)?A t (k 2 ,i)V>(k 2 + q ji ) j (33) 

where g 2 — e 2 /e and q\ = q\ + q 2 . 
Under the scale transformation 



fci,a = 6*1,2, k 3 = b zs k 3 , t=—t, (34) 

b z 

the Coulomb interaction V(q) = g 2 / (q± + q 2 ) transforms as 



-J—d 2 



b 2 g 2 
b g 



if z 3 > 1, 
if z 3 = 1, 



where we display the asymptotic form of the Coulomb inter- 
action in the low energy limit. The resulting scaling relation 
of g 2 is 



(J 2 



g 2 b z - Z3 ifz 3 >l, 
g 2 if z 3 = 1, 

g 2 b z+z 3 -2 ifz3 < 1 



Taking into account of the scaling analysis of So and the rela- 
tions of 23 and z discussed previously, we can see that g 2 is al- 
ways irrelevant except when z = z 3 — 1 where g 2 is marginal. 
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FIG. 6: (Color online) One-loop diagrams that determine the flow of 
the coupling constants. (Left) electron self energy. (Right) photon 
self energy. 



Therefore in the following we perform a one-loop perturbative 
renormalization group analysis to confirm the relevance or ir- 
relevant of the marginal coupling g 2 when z = z 3 = 1. 

To perform the renormalization group analysis we first in- 
troduce a bosonic Hubbard-Stratonovich field ip to treat the 
instantaneous Coulomb interaction. Then the Euclidean ef- 
fective action describing the low energy electrons coupled to 
long range Coulomb interaction can be written as 



SV 



(fx 



^{7o(5o + ig<p) + u(7i9i + 72<9 2 ) - Ad 3 }ijj 



(35) 



where the matrices 70,1.2, satisfying — 2<5y, are de- 

fined as 70 = T3,7i = T2,72 = —T\. tp (ip) represents an 
electron (photon) field operator and ip = ^70- The non- 
interacting electron and photon Green's functions are given 
by 



G (fc ,k) 
D (k ,k) 



~i[7 fc + u(7i fc i + 72^2)] + Ak 2 



kl 



> 2 k\ 



A 2 k\ 



kj + 



(36) 



where k\ = k\ + k\. Here we introduce 77 describing the 
anisotropy of the photon propagators, which is induced by the 
anisotropy of the electron dispersion. (7/ = 1 if there is no 
coupling between the electron and photon.) We perform a 
one-loop Kadanoff- Wilson RG analysis of Share- We trace out 
high energy modes in the energy range of [Ee~ e , E] consider- 
ing constant energy surfaces and obtain a renormalized action. 
Because of the anisotropic dispersion, we have introduced two 
momentum cutoffs A± and A3. 

Let us first consider the renormalization of v and A by com- 
puting the electron self-energy At one-loop order, T,(k) 
can be written as 



where J indicates the integration over the large mo- 

mentum lying within the high energy shell. For convenience, 
we define a dimensionless quantity A = Then S(fc) can 
be written as 

E(fc) =- ivfrih + j 2 k 2 )£h(a,A,r)) ~vA 2 k 3 £I 2 (a, A,r)) 



vA 

a7 



k 2 £I 3 (a,A,r,) 



(38) 



where £ = In b is the usual RG parameter. Ii(a,A,rj) 
(I 3 (a,A,rj)) describes the renormalization of v (A). 
I 2 (a, A, 77) merely describes the shift of the critical pressure 
that is required to tune the system to the quantum critical 
point. Explicitly 11,2.3(0;, A, if) are given by 



h(at, A,r)) = — A 6 / dx 



ol - 3 / 1 — x 



7T J (rjx 2 - A 2 x 4 + A 2 ) 2 



I 2 (a,A,r]) 



a 1 



dx- 



7T V JO (x 



2 _ A 2 x 4 1 A? - 
V ') ' 



tt -r (3zV + A 2 x*r) - x 2 A 2 i 1 ) 
I 3 (a,A, v ) = -A I dx {vx2 _ A2x , + A2)3 ■ 



(39) 



where the fine structure constant a = 4— . 

47TV 

Next we compute the one-loop photon self-energy II(fc), 
which is given by 

n(k)=9 2 J° C fj: J -^Tr[7oGo(g)7oG (fc + <?)], 

(40) 

To determine the renormalization of the bare photon propa- 
gator -Do(fc) by IT(fc), we expand Il(fc) in powers of ki up to 
quadratic order, which leads to 

Il(fc) = -kl£h(a,A,r]) - k 2 lI 5 (a,A,ri), (41) 

in which 

— . 3 a T 8 T 

h(a, A, 77) = -jTp^-, h(a,A,r]) = —aA, (42) 
1U7T A zLir 

Including S(fc) and II(fc), the full effective action can be 
written as 



Srenorm = Stare + / d X1p(-T,)lp + d X-tp(-TL)ip, 



1 



(43) 



The one loop RG equations can be obtained in the following 
-loGo(q)joDo(k — q) way. At first, we rescale the space-time coordinates, 



g 2 f d 3 q [iv( liqi + l2 q 2 ) + Aq 2 ] 



xq = £ob z , x\ = xib, x 2 = x 2 b, x 3 = x 3 b Z3 , 



■2 J {2irf ^/rfqj + A 2 q* 
1 



(44) 



{(k ± -q x ) 2 +r,(ks-q3) 2 Y 



(37) 



where the tilde is used to indicate the scaled variable. Then 
we introduce renormalization constants Z^, Z v , Z g , Z v , Z^, 
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I\ Kia/A for small A, the RG equation for a is given by 



dine 



a 3 a 
4 ~ HhrT 



(47) 



a 



which shows that the velocity renormalization (the first term 
in Eq. d47l >) leads to the logarithmic decrement of a. How- 
ever, due to the exponential decrement of A (it is an irrelevant 
variable), the dielectric constant renormalization (the second 
term in Eq. ( |47| i) dominates the RG flow of a in the small A 
limit, which changes the logarithmic decrement of a to the 
exponential decrement in the low energy limit. Namely, the 
irrelevant variable A plays the role of a dangerously irrele- 
vant variable. In the end, the coupling constant a, which was 
marginal at tree level, turns into an irrelevant variable under 
the one-loop renormalization group correction. 

The singular 1/A behavior of the dielectric constant renor- 
malization I4 can be understood in the following way. When 
A = 0, the photon self-energy is given by 



FIG. 7: (Color online) The RG flow for 2:3 = 1 in (a,A) space. There 
is only one fixed point at the origin. Both a and A exponentially 
decrease to zero. 



— 1/2 7 

and Z v in the following way, ip = Z, ' ip, (p 



,7-1/2 - 
Z v tp, 



_ = Zg 1/2 g,v = Z^d, A = Z A 'A, r) = Z^fj. Then we 
assume that the renormalized action 5 renorm has the same form 
as the bare action Sbare, which leads to 

dhiv 



dl 
din A 

dl 
ding 2 

dl 
din rj 



= z - 1 + Ii(a,A,rj), 
= z- 2z 3 + I 3 (a,A,ri), 
= z - 23 - I 4 (a,A,r)), 



2-2z 3 - I 4 (a, A, rj) + -h(a, A, rj), (45) 
V 



or, in terms of dimensionless variables a = -r— and A = 



, the RG equations are 
din a 



dl _ 
dlnA 

dl 
dim] 



dl 



1 - z 3 - Ix(a, A, rj) - h{a, A, rj), 

1 - 2z 3 + h(a,A,r)) - h(a, A, rj), 

2 - 2z 3 - h(a, A, n) + -I 5 (a, A, 77). (46) 



Now let us describe the RG flow for z 3 = 1 case in which 
the fine structure constant a is marginal. It is interesting to 
note that the RG equations for a consists of two parts, 1\ 
from the velocity renormalization and I4 from the dielectric 
constant renormalization. (Note that since the charge e is pro- 
tected from perturbative renormalization due to gauge invari- 
ance, the renormalization of the coupling constant g 2 = e 2 /e 
results from the dielectric constant renormalization.) Since 



n(fc) = -?-k 2 



9_ 
8v' 



d 3 q 1 



(48) 



Since the integrand is independent of q 3 momentum, the en- 
ergy shell integral simply becomes 



d 3 q = 



which leads to 



dq 3 



d <7j_ = 2 A3 / 2irq±dqi 

*A ± e- e 



16ttV x M 



(49) 



k 2 ± (^)l, (50) 



In the case of the isotropic system such as the 3D Dirac sys- 
tem, -jN- = 1, which leads to the usual logarithmic renor- 



malization of II(fc). However, in the present system, due to 
the strong anisotropy of the electron dispersion, j^- = -i be- 
comes strongly enhanced in the low energy limit, which is the 
origin of the singular enhancement of the dielectric constant 
renormalization. 



The RG flow near the fixed point for 23 = 1 

Since A is an irrelevant variable, it will become very small 
near the fixed point. When A <C 1, the RG equations are 
given by 



dim 



dlnA 



a 3 a 



IOttA' 



dln?7 3 a 8 a A 

dl ~ ~i0^3 + 21tt r) 



(51) 
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5 10 5 10 

I I 

FIG. 8: (Color online) The flow of a (A) varying the initial values of 
A (a) when 23 = 1. (Left) For a fixed initial values of a = 2, the 
initial values of A varies from 0.01 (dotted lines), 0.1, 0.2, 0.5, 1.0, 
2.0, 5.0, 10.0, 20.0, 50.0 (black solid line). (Right) For a fixed initial 
values of A = 1, the initial values of a varies from 0.01 (dotted 
lines), 0.1, 0.2, 0.5, 1.0, 2.0, 5.0, 10.0, 20.0, 50.0 (black solid line). 



The exponential decrement of r\ means that the effective 
Coulomb interaction becomes V(q) « ^- in the end. 

The RG equations for v and A (dimensionfull variables) are 
given by 

—jr- = z - l + Ii(a,A,»7), 

< ^ r = z -2 + h{a,A,rj), (54) 

Notice that the flow equations for v and A depend on the dy- 
namical exponent z. In the low energy limit where z » 1, the 
RG equations become 



For A = A(0)e- e , ■ 



din a 


3 


a 


dl 


" "low 


T 


din 7] 


3 


a 


dl 


' "IOtt 


T 



which leads to 



^ 7i(a,A,r?) = - « const x e \ 
dhxA 

= -1, (55) 



(52) 



which lead to 



= 3 , r. v(e) = v(0)e-™** e - e -+v(Q), A(£) = A(0) e - e . (56) 

A{£) = A{0)e- e 7 

_ Interestingly, the velocity saturates to a constant value while 

— v{fy e > (53) ^4 decreases exponentially. 



